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Abstract 

A random lattice model with dilute interlayer bonds of density p is proposed 
to describe the under doped high-T c cuprates. We show analytically via an 
appropriate perturbation expansion and verify independently by numerical 
scaling of the conductance that for any finite p the states remain extended 
in all directions, despite the presence of interlayer disorder. However, the 
obtained electronic transport is highly anisotropic with violent conductance 
fluctuations occuring in the layering direction, which can be responsible for 
the experimentally observed metallic "in-plane" and semiconducting "out-of- 
plane" resistivity of the cuprates. 



PACS numbers: 72.15.Rn, 71.30.+h, 74.25.Fy 
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In recent years there is a great interest in the problem of Anderson localization for 
anisotropic systems , in connection to the unusual normal-state transport properties of 
the layered high-T c cuprates. The in-plane resistivity for most of these materials exhibits 
metallic behavior, inreasing linearly with temperature over a wide temperature range, while 
the out-of-plane resistivity is reminiscent of semiconductors increasing rapidly at low tem- 
peratures |]7HlO| . This contrasting behavior of the parallel and the perpendicular resistivities 
was observed in Bi2Sr2_ x La a; CuOj / far below T c , down to the lowest experimental temper- 
ature |]TT |. Logarithmic divergencies of the corresponding resistivities accompanied by a 



nearly constant anisotropy ratio are, instead, observed in the underdoped L^-zSr^CuC^ 
suggesting an unusual three-dimensional (3-D) insulator. 

There is currently no consensus on the explanation of the peculiar afr-plane "metallic" and 
the c-axis "semiconducting" transport properties of the cuprates, although it is commonly 
believed that the anisotropic conductivity is, somehow, related to their layered structure. 
It is also widely thought that understanding the peculiar c-axis transport might have im- 
portant consequences for the high-T c superconductivity in the cuprates. In this respect, 
many theoretical models investigate localization in the presence of structural anisotropy, for 
example, the Anderson model with diagonal disorder and anisotropic hopping t between the 
2D planes was recently shown to display very different localization lengths in the parallel 
and the perpendicular directions but the same critical behavior 0. Although this is in 



agreement with the scaling theory of localization for anisotropic systems fl2|--|T5|| is rather 
inconsistent with the observed contrasting resistivities of the cuprates, which has often been 
regarded as evidence for their non-Fermi liquid nature . 



The intense study of the cuprates did not remove many conflicting discrepancies be- 
tween theory and experiment so that a proper description of anisotropic disorder is strongly 
desirable. It must be pointed out that previous disordered site potentials with fixed but 
different hoppings in the lattice axes imply anisotropy only in the band structure and not in 
the disorder configurations which remain isotropic. We propose that in order to understand 
the high-T c cuprates disorder in the anisotropy must be also incorporated. Anisotropic site 



randomness in a form resembling a random superlattice with lateral inhomogeneities gave 
anisotropic localization if the anisotropy was below a critical value, even for arbitrarily small 
disorder || . However, this approach is not a very suitable description for the cuprates either, 
since the CuC>2 planes are believed to be identical with no superlattice-like disorder in the 
layering direction. The important contribution to the c-axis transport in these materials is 
expected to arise from electron scattering in the "insulating" layer between the conducting 
CuC>2 layers 0-0. On the other hand, in almost all high-T c cuprates doping impurities or 
oxygen vacancies always occupy the insulating layers between the conducting "pure" Cu02 
planes which leads to interlayer disorder. 

Having the above considerations in mind we propose a very simple 3D model with 
anisotropy both in the band structure and the disorder configuration. The system con- 
sists of perfect parallel lattice planes randomly connected by interplane bonds described by 
the Hamiltonian 

# = EE |/,m)(/,n|+^t am (|/,m)(/ + l,m|+H. c), (1) 

/ (m,n) i,m 

where m, n are the site indices in every plane and I the plane index. The first term in Eq. (1) 
describes nearest-neighbor intraplane hopping of unit strength which sets the energy unit. 
The second term denotes binary alloy-type of disorder which amounts to interplane hoppings 
t placed at the random planar positions m, implying 100(1 — p)% broken longitudinal bonds 
which inhibit particle migration in the layering direction. The random variable p^ m obeys 
the distribution 

P(p) = (l-p)5{p)+p5(l-p), 

where p is the perpendicular bond density. As a first step towards the explanation of the 
normal-state properties of these materials diagonal disorder is ignored so that the in-plane 
lattices are assumed perfect and the disorder represented by p is due to impurities or oxygen 
vacancies in the insulating layer among the 2D planes. The electrons propagating in this 
system encounter anisotropy in the disorder due to the random interplane links, in addition 
to band structure anisotropy due to t being different than one. 
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A convenient representation for H is the Bloch-Wannier basis 



|k||,/> = -L5>* k n- m |/,m>, (2) 

'-/V|| m 



where kn is the parallel momentum which characterizes the periodic plane states and N\\ the 
number of sites in every plane. In order to investigate the possibility of localization across 
the planes we define the Green function 

G^l-^l'-t) = -i9{t){[c Htl {t)i cl ?f (0)] + ), (3) 

whose diagonal elements give the probability for finding an electron on layer I at time t with 
momentum kn, if initially (t = 0) was on the same layer with the same momentum. The 
Fourier transform of the diagonal Green function reads 

G,(k,,t,k„/;fl) = g _ 1 . (4) 



with the self-energy Sy written as |TT.18| 



mk h l,E) = e ll (k ll ) + Y,Y,Tj n \ (5) 

n=l j 

where the n-th order term is a sum over all paths j starting and ending on the initial layer 
state with 

if = V(k h i; k h ,h) ft V{k E-^) 1+l) ( 6 ) 
and k|| n+1 = k||,/ n+ i = I. The perfect plane Bloch states characterized by k^ have energy 

e[|(k| U ) = 2co8(k Xi ) + 2co8(k Vi ), (7) 

while the interplane hopping which appears in Eq. (6) in the plane-diagonal representation 
takes the form 

V(k lUl k;k lli+1 ,l i+1 ) .i^e^.^l'XMi, (8) 

| m 

with the sum over m denoting only the layer sites having bonds present to the nearest 
neighbor layers. It is seen from Eq. (8) that the diagonal in momentum (ky. = ky ) matrix 



elements are exactly pt and for a given order n if E lies within the 2D band ei^kyj the 
most divergent term of Eq. (5) comes from the path which has intermediate states |kii.,ZA 
with e||(k||.) = E. The corresponding term is (pt/(E — e||(k||J) n and nearest neighbor plane 
states with the same ky are predominantly connected. It can be also seen that momentum 
scattering is always accompanied by interlayer hopping since a change of kn always leads to 
a change of layer. Therefore, H can be tranformed into an "undisturbed" term Hq which 
represents a perfect anisotropic 3D lattice with on-plane (interlayer) hopping 1 (pt), plus 
the rest which forms a "random" Hi hamiltonian part with matrix elements 



v(k,k') = { N EiEm ' 6 -iork^k,,, (g) 

0, for kn = kf„ 



expressed in the i^o-diagonal basis |k) = |k||,fc 2 ), where N is the total number of lattice 
sites, I the layer index, k z the perpendicular momentum and the sum over m/ runs over 
sites which are connected to nearest neighbor planes by interplane bonds. In the adopted 
basis Hi has only off-diagonal matrix elements in the momentum kii between planes since 
the diagonal elements were already included in Hq. 

In order to determine the properties for the eigenstates of H characterized by k, with 
corresponding eigenvalues e(k), we can expand the corresponding self-energy S(k, E) as in 
Eq. (||) via the expansion parameter 

T(E,k') = V (10) 

and if Eq. (9) is substituded into Eq. (10) 

evidently vanishes on average taken over a large number of configurations, since all the on- 
plane site positions will be exhausted for the present interplane bonds by the additional sum 



for the average which leads to cancellations in Eq. ( |TT]) from Bloch theorem. A typical T for 
a given random configuration fluctuates around zero with a variance of um ~ Qcr 2 , where cr 2 is 
the variance of a single term and Q the total number of terms in the sum. According to the 



denominator E — e(k") the most important contribution comes from states with momenta 
on an equal-energy surface in the k-space whose number is of the order of N» = L 2 . The 
total number of present interplane bonds is pL 3 which gives Q ~ pL 5 , so that for a given 
configuration the sum is ~ ^fpVi . Finally, by including the normalization factor iV = L 3 in 
the denominator T(_E,k') ~ ^/pL~^ which guarantees rapid convergence of the self-energy 
in k-space to 

(12) 

The configuration average ( ImJ2k' E-^)-io+ ) * s a ^ so com P u t e d as a function of p and the 
obtained semicircular form can be fitted by p(E)t 2 p(l —p), with E within the H band with 
p(E) the density of states. This allows to estimate the lifetime of states r ~ p ^ t 2 p ^_ p ^ and 
the corresponding mean free paths 

ml 1 

\„ — II ~ i ^ = ^ = ns'i 

11 u p{E)t 2 p{l - p)y/2 + pH' 2 ' y } 

A , = VL ~ P (14) 

by taking into account the Fermi velocites unny It can be seen that for small-p the parallel 
(perpendicular) mean free path is proportional to 1/p (p) so that the interlayer disorder 
affects more the perpendicular transport although cannot localize the states even in this 
direction. 

The predicted extended nature of the wave functions is independently verified by numer- 
ical scaling of the conductance in finite systems. The parallel (perpendicular) conductance 
G\\(±)(L) as a function of the linear size L is obtained for a cubic system at the Fermi energy 



E from the multichannel Landauer-Buttiker formula [IS -21 



G U±) (L) = (e 2 /h)Tr(T+ ±) T u±) ), (15) 

where Tjim is the transmission matrix for electronic propagation along the parallel (per- 
pendicular) direction which can be obtained by transfer matrix techniques. The matrix T± 



turns out to be singular due to the presence of zeros for the missing bonds. This difficulty 
can be dealt with succesfully if we introduce hard-wall boundary conditions in the parallel 
directions and assign a unit incident wave amplitude for one channel, with zero for the rest, 
at a time. The recursion relations for the corresponding wave function coefficients are easily 
established along the parallel direction while along the difficult perpendicular direction the 
assignment made for the incident waves is used as a boundary condition. This is a convenient 
method which allows to consider propagation in the difficult layering direction by a transfer 
matrix product in the easy parallel direction since the L 2 equations for the the unknown 
coefficients can be solved. The absence of missing bonds in the parallel direction makes the 
recursion relations no longer singular and T± can be obtained for the output channels, for 
every assignement made for the input channels. In the calculations we have taken averages 
over up to a 5000 random configurations in each case to suppress fluctuations. 

In Figure 1 we show the scaling behavior of the dimensionless conductance g\\r n(L) = 
G\\(±)(L)/(e 2 /h) with anisotropic interplane coupling t = 0.3 for various bond densities p. 
The parallel conductance increases ballistically (~ L 2 ) for small L and linearly for higher 
L. In the large size (L — > oo) limit the scaling function (3{L) = ding /din L is positive 
for g\\t±), which implies extended states in both directions for any p in agreement with our 
analytical results. However, the transport behavior found is essentially different in the two 
directions. In Fig. 2 we show the energy dependence of the conductance for p = 0.5 and 
t = 0.3 where the ratio of the conductances is close to the estimate g\\/g± — (t\\/t±) 2 for 
t\\ = 1 and t± = pt = 0.15. Another key finding from Fig. 2 is a rather smooth g\\(E) while 
g±{E) displays violent oscillations as a function of E. The dips in g±(E) can be ragarded 
as "minigaps" which might have an effect similar to a semiconductor leading to insulating 
kind of behavior for the out-of-plane conductivity when the Fermmi energy is varied. 

We argue that the proposed model is closely related to high-T c cuprates. As the tem- 
perature increases the inelastic scattering in these materials due to phonons, spin waves or 
other excitations within the CuC>2 layers can cause a decrease of the characteristic inelas- 
tic scattering lengths k n . If the temperature is so high that U n becomes smaller than the 
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mean free path caused from the "random" Hamiltonian H\ the transverse conductivity be- 
comes perfectly metallic. Experiments for Bi 2 Sr 2 CaCu20 8 and underdoped La 2 _ a; Sr :r Cu04, 
YBa2Cu306+x give out-of-plane resistivity which has a semiconductor-like temperature de- 
pendence at low temperatures and a linear-in-T behavior at high temperatures with char- 
acteristic crossover temperature T* between the two regimes, which decreases by increasing 
doping in La^^Sr^CuC^ and YBa 2 Cu306+;c 0-0- We notice that if we relate the bond 
density p with the doping density of the high-T c materials the obtained p-dependence of the 
perpendicular mean free path can be used to explain qualitatively the reported experimental 
behavior. It must be pointed out that the relation between the cuprate doping density and 
the bond density p is very natural, since an increase in the number of the doping impu- 
rities or the oxygen atoms in the layer between two Cu0 2 planes increases the number of 
hopping paths between the two planes. As T* decreases further below the superconducting 
transition temperature T c the out-of-plane normal-state resistivity also becomes metallic, 
which has recently been observed in high-quality single crystals of YBA 2 Cu3C>7 and other 
high-T c cuprates |22[ corresponding to the absence of disorder (p « 1 in the proposed model) 
with almost infinite perpendicular mean free path. This effect can be observed only in the 
perpendicular direction since the parallel mean free path is always much longer and the 
inelasting scattering more dominant in ths direction. 

In summary, we have introduced a simple layered lattice model for high-T c cuprates based 
on the fact that these materials consist of pure Cu0 2 planes while doping impurities or 
oxygen vacancies occupy randomly the interplane hopping paths. We introduce anisotropic 
disorder described by the interplane bond density p, in addition to the usual anisotropic 
band structure expressed via the interplane hopping strength t. In the absence of diagonal 
disorder we show extended states in both directions but the obtained mean free path and the 
conductance in the perpendicular direction is much smaller than the parallel one. Moreover, 
the perpendicular conductance fluctuates strongly as a function of energy, which might 
lead to an insulator-like temperature dependence of the conductivity in this direction. To 
conclude, we have demonstrated anisotropic transport which can account for experimental 

8 



facts of high-T c cuprates. More work is needed in order to understand the corresponding 
metal-insulator transition of layering materials also in the presence of diagonal disorder 
and/or a homogenous magnetic field. 
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FIGURES 



FIG.l. The g\\ as a function of the linear system size L for a cubic system of parallel 
planes with randomly placed interplane bonds of density p and strength t = 0.3. In the inset 
g± is shown for the same system. 

FIG. 2. (a) The energy-dependent g\\ for a layered system with L = 10, 15 and t = 0.3, 
p = 0.5. (b) A much smaller gj_ is shown which exhibits violent oscillations as a function of 
energy. 
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